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1 Introduction 



Recently there has been a surge of interest in lattice actions which are closer 
to the continuum limit than standard lattice actions, in the sense that arti- 
facts due to the finite lattice spacing are suppressed. The hope is that such 
improved actions will allow Monte Carlo simulations to reach the scaling 
region even on rather coarse lattices. 

Based on renormalization group concepts, it has been known for a long 
time that perfect actions, i.e. lattice actions without any cutoff artifacts, do 
exist [|l|]. However, it is very difficult to construct - or even approximate 
- such perfect actions. Recent progress is based on the observation that 
for asymptotically free theories the determination of the fixed point action 
(FPA) is a classical field theory problem ||. At infinite correlation length, 
a FPA is perfect. Away from the critical surface, this is in general not the 
case any more; at finite correlation length, FPAs are referred to as "classically 
perfect actions" . However, they are considered as promising improved actions 
even at moderate correlation length. In particular, the deviations from the 
(quantum) perfect renormalized trajectory are likely to set in only at the two 
loop level H, [5| . Indeed, a drastically improved scaling behavior of the fixed 
point action has been observed in some asymptotically free models || ffl, 
and is especially hoped for in QCD. 

A classically perfect action can be constructed as a fixed point of a block 
variable renormalization group transformation (RGT). Usually one fixes a 
finite blocking factor n and builds block variables on a coarse lattice of spac- 
ing n, determined by a fine lattice with unit spacing. Then one expresses all 
quantities in the new lattice units, i.e. one rescales the coarse lattice back to 
unit spacing. At infinite correlation length and for suitable RGT parameters, 
an infinite number of iterations of this RGT step leads to a finite FPA. For 
most interacting theories, this scheme is the only way to construct FPA's 
non-perturbatively. There one has to perform the blocking transformation 
numerically (and hope for swift convergence). 

For free or perturbatively interacting theories this construction can also 
be achieved in only one step, by a technique that we call "blocking from 
the continuum" . This method proceeds analytically in momentum space. It 
has been applied to construct perturbatively perfect lattice actions for the 
Schwinger model and for QCD, using fixed point fermions of the Wilson type 
[§, [?]]. In this paper we present a number of ingredients for the construction 
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of a classically perfect QCD lattice action with staggered fermions, using the 
same technique. This provides a construction scheme alternative to Ref. [0]. 
A motivation is that staggered fermions are particularly useful for the study 
of chiral symmetry breaking. 

In section 2 we derive the perfect action for free staggered fermions by 
blocking from the continuum. In section 3 we discuss the dispersion relation 
and the pressure of such fixed point fermions with and without truncation. 
In section 4 we show how to block the gauge field consistently and we arrive 
at the corresponding FPA for free gauge fields. This implies a "classically 
perfect" static quark- ant iquark potential, which is compared to the (perfect) 
Coulomb potential in section 5. Section 6 contains our conclusions and an 
outlook on possible applications. 

A synopsis of the results presented here is included in Ref. [§], where we 
also illustrate the perturbative blocking from the continuum and where we 
outline how perturbatively perfect actions can serve as starting points for the 
non-perturbative search for FPA's. 



2 Perfect free staggered fermions 

The fixed point action for free, massless staggered fermions in two dimen- 
sions has been derived before ||, |!|. That derivation used a block factor 3 
RGT, which does not mix the pseudoflavors fl(| , and which could be iterated 
analytically an infinite number of times. We are going to show how one can 
reproduce this result in one step by blocking from the continuum, f] We 
also generalize the result to higher dimensions and to arbitrary masses, so 
we obtain an entire renormalized trajectory for free fermions. 

First we introduce our notation for staggered fermions. We divide the 
lattice into disjoint hypercubes with centers x. Our lattice has spacing 1/2, 
such that these centers have spacing 1. Each hypercube carries a set of 2 d 
fermionic and antifermionic pseudoflavors on its corners, where the dimension 
d is assumed to be even. We denote those Grassmann variables as x p x , xg, 
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This has also been achieved by H. Dilger, in a way which emphasizes the relationship 



to Dirac-Kahler fermions 11 . A similar suggestion occurred earlier in Ref. |12|]. As a 
further approach one could first construct perfect naiv fermions, including all the doublcrs. 
and then build perfect staggered fermions from them. We thank M.-P. Lombardo for this 
remark. 
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where p is the vector pointing from the cell center x to the corner, where the 
variable lives. Hence the standard action for free staggered fermions reads 



Sst[x,x]= E xS[E( r f^ + ^ r 5^)+2m^'] X ^ (2-1) 

where we have defined the following quantities 

T pp' = [§p-M*J + 8P+WJ] (Tll (p), 

r gP' = [5P-A/v_ 5 P+A/2y K(/3)) 

a (x) = { 1 ^»<» x " even 
^ I —1 otherwise 

^P-Xx l^(.Xx+jl Xx— 

ajx£ = x^ +A + xS- A -2xS- (2.2) 

Now we are going to construct a perfect action for free staggered fermions. 
Instead of iterating an RGT with a finite blocking factor, we send the blocking 
factor to infinity and perform only one RGT. This amounts to the "blocking 
from the continuum": one starts from continuum fields and defines lattice 
variables by integrating over the unit hypercube around the corresponding 
lattice site. In particular, for staggered fermions we also have to take care 
of the flavor structure. We have to transform the continuum flavors into 
staggered pseudoflavors and integrate over regions, which depend on the 
pseudoflavor. Those regions are overlapping unit hypercubes. 

It is well-known how to transform fermionic pseudoflavors in the contin- 
uum limit into flavors by a unitary transformation. We assume that the 
inverse of this transformation has been carried out and we start in the con- 
tinuum with space filling fermionic fields ip, ■?/>, composed of pseudoflavors. 
Then we build staggered lattice fermions as 



Xx 



= f d d y V{y) , Xx = f d d y V{y) , (2.3) 



where c z is a unit hypercube with center z. We define 

d 



n(p) = / d*yexp(ipy) = f[ ^ ; p, = 2sin(^/2), (2.4) 

Jc u=l P» 
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and obtain in momentum space 

x P(p) = ip p {p + 27cl)U{p + 27cl)e l{p+27Tl)p 

= E <(p + 27r/)n p (p + 2vr/), (2.5) 
i&z d 

where 

IP(p) -- 

and p is in the Brillouin zone B =] 
choose the following RGT, 

e -s[x,x) _ /" T>ipVipVf]Vr]exp { - s[^,<0] 

+ -L / d d p([x p (-p) ~ E ^("P - 27rl)IF(-p - 27rl)] V P(p) 

+fj»(-p)[ x »(p)- £ n»(p + 27Tl)r(p + 2irl)} 
ie% d 

+fj p (-p)A pp ' V p '{p))}, (2-7) 

where the summation over p, p' is understood. Here s^,^] is the continuum 
action and fj p , rj p are auxiliary staggered Grassmann fields living on the same 
lattice sites as x p an d X p ■ We do not enforce the blocking relation ( |2.5|) by a 
5-function, but only by a smoothly peaked function; the 5-function is smeared 
to a Gaussian by the matrix A. For the latter we make the ansatz, 

A pp' = a5 P P ' _ lC p^T p /e tp{p ' p '\ (2.8) 

where a and c are referred to as mass-like and kinetic smearing parameter. 
They can be tuned to optimize the locality of the perfect action. In particular 
for a = the RGT preserves the remnant chiral symmetry U(l) <8> U(l) (at 
m = 0), which is therefore explicitly present in the perfect action. 
Integrating the RGT ( |2.7| ) we obtain the perfect lattice action 

S[X,X\ = -Jxs \/ d P X p {-p)W(pr l ] pp 'x p '(p), (2.9) 



e ipp U(p) (2.6) 
- 7T, n] d . In analogy to Refs. P, R] we 
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with the free propagator 

A ' ipr ' = ,£ nP(p+M) i fa + 2,th, +m np,(p+2 " )+ - 4W ' 

= -ia„(p)r«" V*"'" + 13^)6'" ', (2.10) 
where a^{j>) and /3(p) are given by 

For p = (pi, 0, . . . , 0), i.e. in the effectively one dimensional case, the sum 
over I collapses to a sum over Ii 6 1, which can be computed analytically. 
In this case, it turns out that for 

sinh(m) — m cosh(m/2) — 1 

a = - ; c = - , [2.12) 

m m z 

the action turns into the standard action, which is ultralocal; the range of 
couplings does not exceed nearest neighbors. In higher dimensions the sum 
over I has to be done numerically and we can only obtain locality in the 
sense of an exponential decay. It turns out that the same choice of smearing 
parameters a and c still yields an extremely local action, i.e. the exponential 
decay is extremely fast, in analogy to our observations for Wilson fermions 
[0. |9| . Hence we are going to use the smearing parameters given in eq. ( f2.12| ) 
in any dimension. 

In the chiral limit m = this result coincides with the FPA obtained in 
Refs. H m . The optimally tuned value of a vanishes in this limit, f\ hence 
we obtain in the massless case an extremely local perfect action, which still 
has the remnant chiral symmetry. 

The inverse propagator can be represented as 

[A^p)- 1 ]^' = p^p)rfe ip(p ~ p/) + X{p)5 pp ' 



3 This is in contrast to Wilson-like fermions, where a chiral symmetry breaking smearing 
parameter is required for a local FPA pi. 
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The function p M (p) is antisymmetric and 2tt antiperiodic in p^. In all other 
momentum components it is symmetric and 2tt periodic, and the same holds 
for X(p) in all components of p. Hence we can expand these functions in 
Fourier series, 

PM = Y,P^ PZ > HP) = E X m ipn , (2.14) 

z ne% d 

where G {±1/2, ±3/2, . . .} and z v 6 2Z, v ^ /i. There is a symmetry 
under permutation and sign flip of A n in all components of n, and of p^ z in 
all z Vy v fi. Furthermore p^ z is antisymmetric in z^. 

Expressed in these quantities, the staggered fermion action takes the form 

s[x,x] = E Z{Ep^( x )l SP ~* /2,p 'xi +z+i > + 8 p+N2 ' p 'xi +z ] 

xe% d p-' z 

+ E A nXS+n}- (2.15) 

n 

The couplings for standard staggered fermions are 

P^z = (#2^,1 - hzp-i) II ^.o ; A n = 2m 5 nfi . (2.16) 



For the perfect action, defined in eqs. (|2.9|) - (|2.12|) , the largest couplings are 
given in Table 1 and 2 for masses m = 0, 1, 2 and 4, and their decay is 
plotted in Fig. 1 and 2. We observe an extreme degree of locality, which is 
very important for practical purposes. For numerical application the perfect 
action has to be truncated to a short range. This truncation ought to alter 
the action as little as possible, in order to preserve the perfect properties to 
a good approximation. 



3 Spectral and thermodynamic properties of 
perfect staggered fermions 

Now we want to address the question, in which sense the action derived in 
section 2 is perfect, i.e. which observables are free of artifacts due to the 
finite lattice spacing. 

First we look at the spectrum. For momentum p = (p,Pa) we see that the 
propagator has a pole at p^ = iy/p 2 ± m 2 , which corresponds to the exact 
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( Za , Zo . Z'i « Za ) 

V I? ^5 05 *± / 


m = 


m = 1 


m = 2 


m = 4 


(1/2,0,0,0) 


0.6617391 


0.5649324 


0.3586038 


0.0730572 


(1/2,0,0,1) 


0.0441181 


0.0335946 


0.0154803 


0.0011324 


(1/2,0,1,1) 

\ / 5 5 5/ 


0.0046569 


0.0032407 


0.0012088 


0.0000711 


(1/2,1,1,1) 

\ / 5 5 5 / 


0.0004839 


0.0003251 


0.0001282 


0.0000114 


(1/2,0,0,2) 

\ / 5 5 5 / 


0.0018423 


0.0012135 


0.0003688 


0.0000063 


(1/2,0,1,2) 

\/ 5^55/ 


0.0001419 


0.0000741 


0.0000096 


-0.0000001 


(1/2,1,1,2) 

\ / 5 5 5 / 


-0.0000264 


-0.0000197 


-0.0000065 


-0.0000001 


(1/2,0,2,2) 

\/ 5^55/ 


-0.0000145 


-0.0000117 


-0.0000046 


-0.0000001 


(1/2,1,2,2) 

\ / 5 5 5 / 


-0.0000126 


-0.0000080 


-0.0000022 


0.0000000 


(1/2,0,0,3) 

V / 5 5 5 / 


0.0000780 


0.0000445 


0.0000090 


0.0000000 


(3/2,0,0,0) 


0.0234887 


0.0172839 


0.0071677 


0.0003303 


(3/2,0,0,1) 


-0.0004933 


-0.0005804 


-0.0004502 


-0.0000389 


(3/2,0,1,1) 


-0.0009913 


-0.0007019 


-0.0002537 


-0.0000067 


(3/2,1,1,1) 


-0.0004819 


-0.0003056 


-0.0000832 


-0.0000009 


(3/2,0,0,2) 


-0.0001210 


-0.0000924 


-0.0000366 


-0.0000008 


(3/2,0,1,2) 


-0.0001011 


-0.0000623 


-0.0000155 


-0.0000001 


(3/2,1,1,2) 


-0.0000462 


-0.0000255 


-0.0000048 


0.0000000 


(3/2,0,2,2) 


-0.0000129 


-0.0000073 


-0.0000015 


0.0000000 


(5/2,0,0,0) 


0.0009439 


0.0006010 


0.0001648 


0.0000018 


(5/2,0,0,1) 


-0.0000600 


-0.0000483 


-0.0000197 


-0.0000003 


(5/2,0,1,1) 


-0.0000577 


-0.0000339 


-0.0000070 


0.0000000 


(5/2,1,1,1) j 


-0.0000187 


-0.0000087 


-0.0000007 


0.0000000 


(7/2,0,0,0) 


0.0000390 


0.0000215 


0.0000039 


0.0000000 



Table 1: The largest kinetic couplings pi(z) for the optimally local, perfect 
staggered fermion of mass 0, 1, 2 and 4. The table contains all the couplings 
> 10~ 5 . 



8 




Figure 1: The decay of the kinetic couplings px(i + 1/2, of the per- 

fect staggered fermion for masses 0, 1, ... 4. We observe that the decay is 
exponential and very fast. For increasing mass it decays even faster. 

continuum spectrum. We note that the II function and the smearing term 
do not affect the singularity structure of the propagator. However, there are 
more poles — one for each / — namely f] 

p it =i\/(p + 27rf) 2 + m 2 . (3.1) 

These poles correspond to higher branches. |] Additional branches are nec- 
essary for perfection, since the lattice imposes 2ti periodicity. 



4 Note that dp4 and the sum over I4 combine to f^L, dp4. 

5 Sometimes higher branches are referred to as "ghosts", which should not be confused, 
however, with Faddeev-Popov ghosts. 
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Figure 2: The decay of the static couplings X(i, in the perfect stag- 

gered fermion action at masses 1, 2, 3 and 4. We confirm that the decay is 
exponential and very fast, especially for large masses. 



As a test case for truncated perfect fermions, it is interesting to observe 
how much harm we do to the spectrum - in particular to its lowest branch 
- if we restrict the couplings to a short range. In Fig. 3 we show the spec- 
trum for a truncated perfect fermion, where we omit all couplings beyond 
±1/2, ±3/2 in the \i direction for p^, and all couplings beyond 0, ±1 in the 
non-p directions of p^ and in all directions of A. The lower branch is real and 
approximates the exact dispersion relation better than the standard stag- 
gered fermion action, and the upper branch is the real part of two complex 
conjugate poles of the propagator. 

An on-shell improved staggered fermion was proposed by S. Naik [|14 
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m = 


m = 1 


m = 2 


m = 4 


(0,0,0,0) 





0.7091500 


0.9799873 


0.5705389 


(0,0,0,1) 





0.0296055 


0.0246899 


0.0034779 


(0,0,1,1) 





-0.0010357 


-0.0015531 


-0.0003674 


(0,1,1,1) 





-0.0012092 


-0.0009711 


-0.0000904 


(1,1,1,1) 





-0.0005181 


-0.0003371 


-0.0000208 


(0,0,0,2) 





0.0009257 


0.0005581 


0.0000181 


(0,0,1,2) 





-0.0000828 


-0.0000703 


-0.0000035 


(0,1,1,2) 





-0.0000567 


-0.0000265 


-0.0000003 


(1,1,1,2) 





-0.0000136 


-0.0000026 


0.0000001 


(0,0,0,3) 





0.0000332 


0.0000133 


0.0000001 



Table 2: The largest static couplings X(z) for the optimally local, perfect 
staggered fermion of mass 1, 2 and 4. The table contains all the couplings 
> 10" 5 . 



By changing the kinetic term of the standard action ([2.1 6|) to 

/9 1 \ 

/W = (g [<Wi - S 2Zfl -i] - ^t[^ m ,3 - II , ( 3 - 2 ) 

he removed the 0(a 2 ) artifacts for the free fermion. These couplings have no 
similarity to our perfect fermion. For comparison, the resulting dispersion 
relation is also shown in Fig. 3. By construction, this dispersion relation is 
good at \p\ << 1, but at \p^\ ~ 1 the lower branch is hit by an upper branch, 
and then they turn into two complex conjugate poles, the real part of which 
is shown in the figure. 

With respect to certain other quantities, even our untruncated action is 
not exactly perfect. When we performed the Gaussian integrals in the RGT 
( |2.7| ), we did not keep track of "constant factors" , which do not depend on the 
lattice fields. However, such factors may depend on other quantities such as 
the temperature. If those quantities are important - which was not the case 
in the spectrum - we can not expect perfect observables. Then one notices 
that the lattice fields are not completely renormalized. The possibility to 
keep track of such factors in a block spin RGT was explored for the Ising 
model in Ref. |I5[ . 
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perfect 
Naik 

truncated perfect 

standard 



1.0 2.0 3.0 

p along (1,1,0) 



Figure 3: Dispersion relation for truncated perfect, Naik and standard, stag- 
gered fermions along (1,1,0). 

As an example we consider the pressure of free, massless fermions in 
infinite volume. In the continuum the relation 



P 

y4 



III 

180 



0.3838 



(3.3) 



is known, where P is the pressure and T the temperature (Stefan-Boltzmann 
law). In Fig. 4 we plot the ratio P/T 4 for the staggered FPA and for the 
Naik and standard and staggered action at N t = 2, 4, 6 ... 40 lattice points in 
the 4-direction (N t must be even to accommodate a complete set of pseud- 
oflavors). The standard staggered action scales better than Wilson fermions, 
in agreement with the fact that their artifacts are of the second order in the 
lattice spacing, whereas the Wilson action is plagued by linear artifacts 

6 For Wilson fermions the ratio P/T A has a peak ~ 1.9 at N t ~ 2, 3 ||. For a discussion 
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We also see that even for the FPA this scaling quantity differs from the 
exact continuum result. This deviation is due to a missing temperature 
dependent renormalization factor. However, the figure shows that this fac- 
tor is typically very close to 1. It differs significantly only on immensely 
coarse lattices. Hence the scaling is considerably improved with respect to 
the standard lattice actions. This is still true after truncation (the same we 
used when discussing the spectrum), although the behavior gets somewhat 
worse. Truncation causes an overshoot, which indicates that the fermion 
moves closer to the standard formulation. For the Naik fermion we confirm 
some improvement as well. Its thermodynamic behavior has also been dis- 
cussed in Ref. [pjj. The dip in the beginning is very similar to the behavior 
of the D234 action [T8|] (the corresponding thermodynamic plot is given in 
|| ) . Note that also the construction of those two on-shell improved fermions 
- Naik for the staggered and D234 for the Wilson type fermions - is very 
similar: in both cases additional couplings are added on the axes (which is 
not too promising for the restoration of rotational invariance). Finally, the 
spectrum is very similar too, the lowest branches are hit by a doubler and 
turn complex around [p M | ~ 1. 



4 A free fixed point gauge action consistent 
with staggered fermions 

Perfect staggered fermions can only be used in QCD if we are able to couple 
them to fixed point gauge fields in a perfect way. This requires a blocking 
scheme for the gauge field, which is consistent with the blocking for staggered 
fermions. The consistency condition is based on gauge covariance. 

Let's go back to a finite blocking factor n. When we block from a fine 
to a coarse lattice, it is convenient to fix a gauge for the fine lattice fields. 
However, this gauge fixing should be restricted to one block, in order to avoid 
long distance gauge dependence. 

Following Ref. |10| we block such that all contributions to a coarse lattice 



variable have the same pseudoflavor, as it is illustrated in Fig. 5. If a coarse 
pseudoflavor lives on a site x', then we consider the hypercube n d with center 



of the artifacts in standard staggered fermions and an improvement program for its matrix 
elements, see [16|. 
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Figure 4: The thermodynamic scaling quantity P/T 4 for free massless stag- 
gered fermions with the standard action, the fixed point actions with and 
without truncation, and the Naik action. 

x', and all the n d fine variables of the same pseudoflavor contribute to that 
block variable. Thus each fine variable contributes to one coarse variable 
of the same pseudoflavor. (In the previous section we considered the limit 
n — > oo of this RGT.) 

A coarse gauge field A! , in terms of fine fields A^ x , which is consistent 
with gauge covariance, is 

b n 1 n_1 

A 'n,x'+n(l/2 = — A - A V,x+(2j+l)p,/2, (4.1) 
xGx' j=0 

where b n is a constant renormalization factor. Its value is chosen such that 
we obtain a finite FPA for the free gauge field. Essentially, b n neutralizes 
the constant factor from rescaling A' , hence dimensional reasons suggest 
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Here our fine lattice has spacing 1/2 and the coarse one n/2, (n odd). 
The sum x G x' runs over all n d fine lattice points that contribute to the 
coarse pseudoflavor living at the site x', and the lattice gauge fields live on 
the link centers. We sum over the straight connections of corresponding fine 
lattice pseudoflavors contributing to adjacent coarse lattice variables. This 
construction is illustrated for n = 3 in Fig. 5. 
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Figure 5: The staggered blocking of fermions and gauge fields, which main- 
tains the pseudoflavor structure and gauge covariance. This figure illustrates 
the case d = 2, n — 3. Each symbol represents a pseudoflavor. Block vari- 
ables are encircled and pick up contributions of the fine variables (with the 
same flavor) in the dashed box around them. The block link (shaded) is built 
from all the links covered by arrows. 

In momentum space this condition reads 



n 



m = KM. n M ( ) = 4sin(p,/4) 



(4.2) 



where the summation extends over V G {0, 1, 2, . . . , n - l} d . Note that the 
momenta of the fine fields are in the zone B 4n =] — 2ir, 2ir] d , and those of the 
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coarse fields in B^/ n =] —2n/n, 2-n /n] d . Since we defined the gauge variables 
on the link centers, A^(p) is 4n antiperiodic in p^. 

Also here we want to send the blocking factor to infinity and block from 
the continuum. We start from a space-filling set of pseudoflavors, together 
with a continuum gauge field a^. In the limit n — > oo the above condition 
(|4. 1|) turns into a Riemann integral of the form 

J d d yM^y)a^x-y) 

r M{ Vll ) \y v \ < 1/2, u^n 
I otherwise 

f i \y,\ < i/4 

3/2-2|^| 1/4 < |^| < 3/4 . (4.3) 
[ otherwise 

This convolution is illustrated in Fig. 6. It is manifestly gauge covariant. Due 
to its architecturally interesting shape we call the function M M the mansard 
function. 




Figure 6: The "mansard blocking" of a gauge field from the continuum. 
Note that 

J d d yM,(y) eMipy) = ITf (p), (4.4) 
and therefore condition (|4.2| ) turns into 

MP) = E a M (p + 47rZ)ITf (p + 47rZ)(-l)^. (4.5) 

The following property can be understood from the above derivation, which 
started from a discrete RGT: if we perform a gauge transformation on an 



MM = 
M(y») = 
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Abelian gauge field, a^(y) — > a^iv) + d^(p(y), then the lattice gauge field 
transforms as A^ x -> A^ x + $ s _ A / 2 - ^ x +(i/2, where $ x = d d xip(x) (c x 
being a unit hypercube with center x). This shows gauge covariance and 
consistency with the fermionic blocking. The same properties were achieved 
in Refs. || [7] for fermions of the Wilson type, where the blocking was done 
by stepwise integration over disjoint cells and the gauge field was convoluted 
with a roof shaped function instead of the mansard function. 

Now we construct the FPA, for free gauge fields , which is adequate 

for staggered fermions. We impose the Landau gauge on the continuum gauge 
fields and choose the RGT 

exp{-SL4J} = J VaVDexp { - J d d p^(-p)p 2 a^p)} 

\ r \ 2 

x exp{-^p d d p-Zy-p)Kp)+ 7 (p)(p M /2) ]D M (p) 
+j£>„(-p)[A /1 (p)- a M (p + 47r/)nf(p + 4vr/)(-l)^]},(4.6) 

where the measure Da contains the gauge fixing factor 5{J2^p^a^{p)) Again 
we impose the blocking condition ( [4.5| ) only by a Gaussian involving a mass- 
like and a kinetic smearing term. is an auxiliary lattice field, defined 
on the same links as A^, in analogy to the fermionic auxiliary fields 77,77 in 
eq. Performing the RGT O yields 

A » = E , \ ,, 2 n; f (p + 4< + a (p)+ 7 (p)(p72) 2 . (4.7) 

This is the FPA in a special gauge that we call "fixed point lattice Landau 
gauge", where A^ obeys 

Esin^( P )-M At (p) = 0. (4.8) 

We follow the standard procedure to rewrite the FPA in a gauge invariant 
form and arrive at 

s[A) = j^J B J d p\M-p)^Apy l Mv), 
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sin (p u I A) A?, ( p) ~ 1 A£ (») ~ 1 sin /4) 
£ A sm (p A /4)A^(p) 1 

This result is similar to the one obtained for the roof- like blocking in |6|, PR. 
There it was possible to choose the smearing parameters such that the FPA in 
d = 2 turned into the ultralocal plaquette action. This can also be achieved 



for the "staggered FPA" given in eq. ( 4.9 



In d = 2, the standard lattice Landau gauge implies 

_ sm( Pu /4)F(p) 
A>)-^ 4s . nW4) , (4.10) 

where F is the plaquette variable defined on the plaquette centers. Inserting 
this into eq. (|4.9|) for d = 2, we find 



^ = 7TV2 f d 2 p]-F(-p)p(p)F(p), 
{2-nY Jba„ 2 



(2tt) 2 Jb 47I 

pip)- 1 = 16[sin 2 |Af(p) + sin 2 |A?(p)] 

= cos^cos^(l-i(p77 2 ) 2 ))(l-i(p72) 2 ; 

+Aa{p)[{p7/2) + (p72) 2 ] + 8 7 (p)(p72) 2 (P272) 2 - (4.11) 



If we choose 

/ x 1 1 2 Pi 2P2 

a (p) = cos — cos — 

w 16 24 4 4 

/ \ 1 1 2 Pi 2 P2 1Q x 

TIP) = cos — cos — , 4.12 

rw 128 288 4 4 ' v ; 

then we obtain p _1 (p) = 1, as desired. Taking this as a guide for higher 
dimensions suggests that the optimal choice for the RGT parameters is 

- 2 

It is important that a(p) + 7(p)(p At /2) is always positive. This ensures that 
the functional integrals in the RGT (|4.6|) are well defined. 
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Indeed it turns out that this RGT yields a very local FPA in d — 4. The 
largest couplings are given in table 3 and the exponential decay is plotted in 
Fig. 7. 




Figure 7: The exponential decay of the free gluon FPA for a RGT of the 
mansard type. 

The dispersion relation of the transverse mansard gluon is again perfect, 
as it is the case for the perfect free fermion. If we truncate the couplings 
to distances < 3/2, we obtain the dispersion shown in Fig. 8. Due to our 
optimization of locality, the spectrum of the truncated fixed point gluon 
approximates the continuum spectrum still very well for small and moderate 
momenta. In practice, since these gauge fields are designed to couple to 
staggered fermions, the relevant momentum region is \p^\ < ir, where the 
dispersion is quite good. We will find further evidence for this from the 
calculation of the quark- ant iquark potential. 
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( Za , Zo . Zi* Za ) 

V 17 ^5 07 *± / 


Pn (z) 




( Za . Zo . Zi . Za ) 

V 1 5 Z 7 O 7 *± / 


Pao(z) 
r lz \ / 


(0,0,0,0) 


3.8112883 




(1/2,1/2,0,0) 

V / 7 / 7 7/ 


0.6609129 


(0,0,0,1) 


-0.3402240 




(1/2,1/2,0,1) 

\ / 7 1 7 7/ 


0.0742258 


(1,0,0,0) 


-0.1541891 




(1/2,3/2,0,0) 

\ / 7 / 7 7/ 


-0.0271129 


(0,0,1,1) 


-0.1150765 




(1/2,1/2,1,1) 

V / 7 / 5 7/ 


0.0148817 


(1,0,0,1) 

\ 7 7 7 / 


0.0508985 




(1/2,3/2,0,1) 

\ / 7^/ 7 7/ 


0.0047334 


(0,1,1,1) 


-0.0476001 




(3/2,3/2,0,0) 


0.0088456 


(1,0,1,1) 


-0.0058551 




(1/2,3/2,1,1) 

V / 7 / 7 7 / 


0.0015899 


(1,1,1,1) 

\ 7 5 7 / 


-0.0077244 




(3/2,3/2,0,1) 

V / 7 / 7 7/ 


-0.0003598 


(0,0,0,2) 


0.0054232 




(3/2,3/2,1,1) 

\ / 7 / 7 7/ 


-0.0010820 


2,0,0,0 


0.0084881 




(1/2,1/2,0,2) 


-0.0041648 


(0,0,1,2) 


-0.0012443 




(1/2,5/2,0,0) 


0.0015083 


(1,0,0,2) 


-0.0065388 




(1/2,1/2,1,2) 


-0.0001960 


(2,0,0,1) 


-0.0039935 




(1/2,3/2,0,2) 


0.0004785 


(0,1,1,2) 


-0.0006582 




(1/2,5/2,0,1) 


-0.0005171 


(2,0,1,1) 


0.0011516 




(3/2,5/2,0,0) 


-0.0007231 


(1,1,1,2) 


0.0004326 




(3/2,3/2,0,2) 


-0.0002568 


(2,0,0,2) 


0.0011159 




(3/2,5/2,0,1) 


0.0001332 


(1,0,2,2) 


0.0005143 




(3/2,3/2,1,2) 


-0.0002152 


(1,1,2,2) 


0.0003562 




(1/2,1/2,2,2) 


0.0001109 


(1,2,1,1) 


0.0002912 




(5/2,5/2,0,0) 


0.0001907 


(0,0,0,3) 


-0.0005163 




(1/2,1/2,0,3) 


0.0002217 


(3,0,0,0) 


-0.0005617 




(1/2,7/2,0,0) 


0.0001005 


(3,0,0,1) 


0.0003064 




(1/2,3/2,0,3) 


-0.0001036 



Table 3: The largest couplings for the fixed point action of the free gluon 
with respect to the "mansard RGT". The table includes all couplings in pu 
with values > 2.5 • 10~ 4 and all couplings in p i2 with values > 10~ 4 . 
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10.0 



fixed point 

Truncated Fixed Point 




Figure 8: The dispersion relation of the truncated mansard fixed point gluon. 

5 Polyakov loop and the static quark- ant iquark 
potential 

It is also possible to construct perfect operators for a given RGT. One adds a 
source term to the continuum action (or the fine lattice action) and includes 
this term perturbatively in the blocking process. In Ref. the perfect XX 
operator was constructed in this way. In principle this method allows us 
to build any perfect composite operators, but in practice this tends to be 
difficult. 

It is easier to construct "classically perfect" fields and operators. Classi- 
cally perfect fields are obtained by minimizing the continuum action together 
with the blocking transformation term. If one inserts h in the RGT expres- 
sion (|4.6|) , then one obtains them in the limit H — > 0. 

For the free gauge field — and the RGT considered above — the classically 
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perfect field reads 



a» = inf (p)A«(p)-%(p), (5.1) 



where is in the Landau gauge and in the fixed point lattice Landau 
gauge. Note that classically perfect fields are defined in the continuum. 

We obtain "classically perfect operators" from continuum operators if 
we replace the continuum fields by classically perfect fields, which are then 
expressed in terms of lattice fields. For instance, we can build a classically 
perfect Polyakov loop as 

(j) c (x) = / dx d a d (x,x d ). (5.2) 



Inserting the gauge field given in eq. (|5.1| ), we obtain in momentum space 

np) = o)A§(a oy 1 am o). (5.3) 

Also this operator is defined in the continuum, in contrast to the standard 
lattice Polyakov loop = A d (p, 0). Of course we can restrict the classi- 
cally perfect Polyakov loop to lattice points by imposing 2n periodicity, 

nP)= E ^ \ r nf(p + 27rr,0)A%p,oy 1 A d (p,0). (5.4) 

In d = 2 we found A|(pi, 0) _1 = p\ , which leads to 

K = \{**i+2** + (5.5) 

i.e. the classically perfect Polyakov loop is ultralocal on the 2d lattice. This 
confirms that the RGT parameters introduced in Sec. 4 optimize locality. 
The correlation function of two classically perfect Polyakov loops, 

yields the static quark- ant iquark potential 

V(#0 = " (2^5=1 / dd ~ l Pjhf^^ °) aA S® ^P(¥r). (5-7) 
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In d = 2 this integral diverges as it stands. We subtract an infinite constant 
such that we obtain the correct behavior at large r. Then integration by 
contour techniques leads to 

r/2 r > 1.5 

(-32r 5 + 240r 4 - 720r 3 + 108O 2 - 330r + 243) /960 1 < r < 1.5 
(32r 5 - 80r 4 - 80r 3 + 440r 2 - lOr + 179) /960 0.5 < r < 1 
(64r 5 - 160r 4 + 400r 2 + 178) /960 < r < 0.5 

(5.8) 

This potential is shown in Fig. 9. It is very smooth (four times continuously 
differentiable) and coincides with the exact potential at r > 1.5. 




0.0 !«- 1 1 1 1 1 1 1 

0.0 0.5 1.0 1.5 2.0 



Figure 9: The classically perfect static quark- ant iquark potential in d = 2. 

In d = 4 the Fourier transform (|5.7| ) has to be calculated numerically. 
The result is shown in Fig. 10 and compared to the static potential based 
on the standard lattice Polyakov loop, which is only defined if r is a lattice 
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vector. For the classically perfect potential, we observe a faster convergence 
of the scaling quantity |f|V(f) to the exact value —1/Att. Also rotational 
invariance is approximated much better, even down to \r\ < 1. The reason 
for the remaining artifacts in the classically perfect potential is related to the 
reason for imperfectness of the pressure discussed in Sec. 3. Here we also 
ignore "constant factors" in the Gaussian integrals of the RGT by taking 
just the minimum of the exponent. The remaining artifacts are exponentially 
suppressed. 



0.10 




Figure 10: The classically perfect static quark-antiquark potential in d — 4 
on an axis and on the space diagonal, compared to the potential obtained 
from Wilson's standard plaquette action. 
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6 Conclusions and outlook 



We have constructed perfect free staggered fermions and perfect free gauge 
fields, which can be coupled consistently. This forms a basis for systematic 
perfect lattice perturbation theory, leading to a perturbatively perfect lattice 
action for QCD with staggered fermions. The next step is the evaluation of 
the perfect quark-gluon and the three gluon vertex functions as it was done 
for Wilson-like fermions J7|, || 19 1. In momentum space these vertex functions 



have almost the same form as written down in J7[ , with some modifications for 
staggered fermions and mansard gluons, which are obvious from the present 
paper. The lattice action thus obtained would be perfect for weak couplings. 
It may already be useful in lattice simulations even at moderate correlation 
length, since the classically perfect action is also expected to be one-loop 
perfect. The next step is to test the actions in simple situations. For example, 
one could study heavy quark physics, even though staggered fermions are 
not designed for such applications. In such a study, the worst artifacts of the 
standard lattice formulations originate from the large mass of the quarks. 
However, they are considerably suppressed for actions which are perfect at 
weak couplings. For Wilson fermions, we observed that even the use of 
the free perfect fermion action, together with standard lattice gauge fields, 
improves the mesonic dispersion dramatically || — although a large additive 
mass renormalization occurs. For staggered fermions, only a multiplicative 
mass renormalization will occur due to chiral symmetry. 

An important aspect of staggered fermions is that they break the flavor 
symmetry when interacting with gauge fields. In lattice QCD with staggered 
fermions there is only one true "Goldstone pion". The other pions remain 
massive in the chiral limit, which manifests the flavor symmetry breaking 
f2~0|j . This is a major obstacle in studying chiral symmetry breaking at finite 
temperatures, since the number of flavors plays an important role in such 
studies. In an experiment with a "fat link" - consisting of a link plus 
staples with varying staple weight — it was possible to reduce the mass of 
the remaining pions by about a factor two |nj], which shows that there is 



a large potential for improvement by using non standard actions. However, 
the use of the Naik fermion does not seem to help here [P^j . When we 
couple perfect staggered fermions and gauge fields consistently, we expect 
the flavor symmetry breaking to be strongly reduced. The "fat link" would 
naturally be incorporated. This could help to make the studies of chiral phase 
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transitions at finite temperature more realistic. As a simpler experiment 
one could combine the truncated perfect staggered fermion with a standard 
lattice gauge field. Perhaps this already helps to decrease the flavor symmetry 
breaking. 

In principle, if the dynamics of the problem involve large gauge fields, 
as is obviously true for the physics of the light pions, it is unclear if the 
improvements obtained perturbatively are sufficient at the typical couplings 
where simulations can be performed. The couplings obtained perturbatively 
could undergo further renormalizations and new couplings may arise.[] In 
such a situation one has to resort to non-perturbative techniques in order to 
obtain the classically perfect action. In this case one rescales the QCD action 
with the gauge coupling g and looks for a fixed point - i.e. a classically perfect 
action - of the rescaled theory. This can be done numerically by real space 
RGT steps, using the parameters for RGTs with finite blocking factors given 
in the appendix. The perturbatively perfect action may serve as a promising 
point of departure for this iteration. Thanks to the asymptotic freedom at 
g = 0, such RGT steps only consist of minimization of the action on the fine 
lattice together with the transformation term. No (numerical) functional 
integration needs to be performed, which simplifies the task enormously. In 
practice one starts from an updated configuration on the coarse lattice and 
determines the minimizing fields on the fine lattice. In this way, the effects 
of strongly fluctuating field configurations can be incorporated. 
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7 However, we would like to emphasize that the perturbatively perfect action will ap- 
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Appendix 



A Parameters for RGTs with finite blocking 
factors and the same fixed point 

In the limit g = (where we determine the FPA), the quarks decouple from 
the gauge fields, hence we can start by constructing a pure gauge FPA. Here 
we provide the analytic ingredients for this purpose. 

In Sec. 4 we derived a FPA for free gauge field, which is consistent with the 
gauge requirements of staggered fermions. There we blocked from the contin- 
uum, but of course the same FPA can also be obtained using finite blocking 
factor RGTs. The optimal RGT parameters, which provide ultralocality in 
d = 2 and extreme locality in higher dimensions, depend on this blocking 
factor n (n odd). So far we only gave those parameters at n — > oo. In this 
appendix we are going to identify them for general n. This is needed for the 
nonperturbative search of fixed points of non-Abelian gauge fields. There one 
does RGTs numerically, and is therefore restricted to finite blocking factors. 
In practice one would choose the smallest value n — 3. 

Assume we are at the fixed point of Sec. 4 in the fixed point lattice Landau 
gauge, and we perform a block factor n RGT, 



where V G {0,1,2, ... ,n — l} d as in condition ( |4.2|) , which is (smoothly) 
implemented here. For the smearing term we need an ansatz, which 
involves more parameters than it was the case for the blocking from the 
continuum, 



exp{-S'[A'}} = VAVDeMS[A,}} 




+*D,(-p) [^(p) MP + ^l'/n)K(p + 4*l'/n){-l)l] }, (A.l) 



ty*,n(p) 



a n + 4sin 2 (np At /4)7 n + sin 2 (np fl /2)uj. 



d 



+ [5 n + 4sin 2 (np M /4)(T n ] J[ cos 2 [np v j 4). 



(A.2) 



27 



The RGT parameters b n , a n , 7 n , u n , 5 n and a n can now be determined 
from the condition that the action be invariant under this RGT. Doing this 
integral and rescaling the coarse lattice momenta into the full zone -84^, we 
obtain 

A9 » = 4E^((P + 47r O/^)n^((p + 4vr/')/n) 2 + ^, n (p).(A.3) 
n v 

We combine the sum over Airl'/n with the sum over 4nl, I e Z5 d , which is 
intrinsic in A 9 , to a sum over 4n£/n, i G ^ a! . After doing some lengthy 
algebra and evaluating a number of trigonometric sums, we find that the 
fixed point condition A^'(p) = A 9 (j)) is fulfilled if 



K = n 



d-2 

n 2 -l 



in 2 



(n 2 - l) 2 

^ = 9 ^ = ■ < A - 4) 

We see that Q^ n (p) is always positive, such that JVD is well-defined. In 
particular, we confirm the result for b n , which we anticipated in Sec. 4 by 
dimensional reasons, and in the limit n — > 00 we reproduce the parameters 
for the blocking from the continuum obtained in eq. (|4.12|). As a further 



check we notice that n — 1 yields a trivial identity transformation, as it 
should. 
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